
IvetenskapB KUNGL 

' °™ i Tekniska 



HOGSKOIAN 



Nonlinear Optics 5A5513 (2003) 
Lecture notes 



Lecture 2 



Nonlinear polarization density 

From the introductory perturbation analysis of the aU-classical anharmonic oscillator in the pre- 
vious lecture, we now a priori know that it is possible to express the electric polarization density 
as a power series in the electric field of the optical wave. 

Loosely formulated, the electric polarization density in complex notation can be taken as the 
series 

-P/i(Wa) = eQ[x^ill{-UJcr\UJa)Ea(uJ„) + x|fj^ ("W^ ; Wi , W2)-Ea (wi (a;2) 

^ ' ^ ^ ' 

^ ^ ^ 

where we adopted Einstein's convention of summation for terms with repeated subscripts. (A more 
formal formulation of the polarization density will be described later.) 

Symmetries in nonlinear optics 

There are essentially four classes of symmetries that we will encounter in this course: 
[1] Intrinsic permutation symmetry: 

x|fa/3-y(-W<T;Wi,W2,C^3) = X^f/3a^(-W<T;W2,Wl,W3) 
= x[37"(-'^'^;^2,W3,Wi) 

that is to say, invariance under the n! possible permutations of {a}^,uj}^), fc = 1,2, . . . ,n. This 
principle applies generally to resonant as well as nonresonant media. {Described in this lecture.) 
[2] Overall permutation symmetry: 

x2/37(-^a; Wi, W2, 1^3) = x2/3-y(wi; -Wa, UJ2,U)3) 
= Xa/3^7K;^2,-W^,W3) 

= xi%^^K:'^2.'^3,-t^a), 

that is to say, invariance under the (n + 1)! possible permutations of (/x, — Wq-), (a^, w^), k = 
1, 2, . . . , n. This principle applies to nonresonant media, where all optical frequencies appearing 
in the formula for the susceptibility are removed far from the transition frequencies of the 
medium. {Described in lecture 4.) 
[3] Kleinman symmetry: 

x|fi/37(-^<^;'^l'^2,W3) = xflp-y{-0Ja;0Ji,UJ2,U>3) 
= xiX7(~'^<^''^l''^2,W3) 
= Xaj-y^(-t^a;Wi,W2,W3), 
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that is to say, invariancc under the (n+1)! possible permutations of the subscripts /i, ai, . . . , a„. 
This principle is a consequence of the overall permutation symmetry, and applies in the low- 
frequency limit of nonresonant media. 
[4] Spatial symmetries, given by the point symmetry class of the medium. {Described in lecture 6.) 

Conditions for observing nonlinear optical interactions 

Loosely formulated, a nonlinear response between light and matter depends on one of two key 
indgredients: either there is a resonance between the light wave and some natural oscillation mode 
of the medium, or the light is sufficiently intense. Direct resonance can occur in isolated intervals 
of the electromagnetic spectrum at 

• ultraviolet and visible frequencies (10^^ s~^) where the oscillator corresponds to an electronic 
transition of the medium, 

• infrared (10^"^ where the medium has vibrational modes, and 

• the far infrared-microwave range (10^"'^ s^^), where there arc rotational modes. These interac- 
tions are also called one-photon processes, and are schematically illustrated in Fig. 3. 

\b) 



1«) 

The one-photon transition. 
Figure 3. Transition scheme of the one-photon process. 

The lower frequency modes can be excited at optical frequencies (lO'^^ through indirect reso- 
nant processes in which the difference in frequencies and wave vectors of two light waves, called the 
pump and Stokes wave, respectively, matches the frequency and wave vector of one of these lower 
frequency modes. These three-frequency interactions are sometime called two-photon processes. 
In the case where there "lower frequency" mode is an electronic transition or in the vibrational 
range (in which case the Stokes frequency can be of the same order of magnitude as that of the 
light wave), this process is called Raman scattering. 
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The Stokes Raman transition. The Anti-Stokcs Raman transition. 

Figure 4. Transition schemes of stimulated Raman Scattering. 

The stimulated Raman scattering is essentially a two-photon process in which one photon at uii 
is absorbed and one photon at uj2 is emitted, while the material makes a transition from the 
initial state \a) to the final state |c), as shown in Fig. 4. Energy conservation requires the Raman 
resonance frequency (electronic or vibrational) to satisfy h^ca = h{u}i — W2), and hence we may 
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classify the Stokes and Anti-Stokes transitions as 

TiD,ca > <^ Stokes Raman 

h^lca < <^ Anti — Stokes Raman 

When the lower frequency mode instead is an acoustic mode of the material, the process is instead 
called Brillouin scattering. 

As will be shown explicitly later on in the coiu'se, optical resonance with transitions of the 
material is an important tool for "boosting up" the nonlinearities, with enhanced possibilities of 
applications. However, at single-photon resonance we have a strong absorption at the frequency 
of the optical field, and in most cases this is a non-dcsirablc effect, since it decreases the optical 
intensity, even though it meanwhile also enhances the nonlinearity of the material. 

However, by instead exploiting the two-photon resonance, the desired nonlinearity can be sig- 
nificantly enhanced whilst at the same time the competing absorption process can be minimized 
by avoiding coincidences between optical frequencies and single-photon resonances. 

One general drawback with the resonant enhancement is that the response time of the polar- 
ization density of the material is slowed down, affecting applications such as optical switching or 
modulation, where speed is of importance. 

Phenomenological description of the susceptibility tensors 

Before entering the full quantum-mechanical formalism, we will assume the medium to possess a 
temporal response described by time response functions R{t) . 

The very first step in the analysis is to express the electric polarization density, which in classical 
electrodynamical terms is expressed as the sum over all M electric charges qk in a small volume V 
centered at r, 

1 ^ 
fe=i 

as a series expansion 

P(r, t) = P W (r, t) + P(i) (r, t) + P^^) (r, i) + . . . + p(") (r, t) + . . . , 

where P^^)(r, i) is linear in the electric field, P^^)(r, i) is quadratic in the electric field, etc. The 
field-independent P(°)(r, t) corresponds to eventually appearing static polarization of the medium. 
It should be emphasized that any of these terms may be linear as well as nonlinear functions of 
the electric field of the optical wave; i. e. an externally applied static electric field may together 
with the electric field of the light interact through, for example, P^^^{r,t) to induce an electric 
polarization that is linear in the optical field. 



R{t - t) 



R = 



t = T t 

Figure 9. Schematic form of a possible response function in time domain. 

In the analysis that now is to follow, we will a priori assume that it is possible to express the 
polarization density as a perturbation series. Later on, we will show how each of the terms can be 
derived from a more stringent basis, where it also will be stated which conditions that must hold 
in order to apply perturbation theory. 
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Linear polarization response function 

Since p(^)(r, t) is taken as linear in the electric field, we may express the linear polarization density 
of the medium as being related to the optical field as 

/oc 
rW(t;r)£;,(r,T)dr, (7) 
-oo 

where r^aHt;r) is a rank-two tensor that weights all contributions in time from the electric field 
of the light. A few required properties of T^a (i; t) can immediately be stated: 

• Causality. Wc require that no optically induced contribution can occur before the field is 
applied, i. e. Tj^2 {t; r) = for i < r. 

• Time invariance. Under most circumstances, we may in addition assume that the material 
parameters are constant in time, such that P^^^(r,t') is identical to the polarisation as induced 
by the time-displaced electric field Ea{r,t'). 

The second of these properties is essentially a manifestation of an adiabatically following change 
of the carrier wave of the optical field. 

By using Eq. (7), the time invariance of the constitutive relation gives that 

/oo 
7iL'(i + ^o;T)S„(r,T) dr 
-OO 

= {Should equal to polarization induced by ^^(r, r + to)} 

/oo 
Tj,^J{t;T)E^{r,T + to)dT 
-OO 

= {r' = T + io} 

/OO 
T^]I{t;r' -to)E^ir,T')dT', 
-oo 

from which we, by changing the "dummy" variable of integration back to r, obtain the relation 

Tl,]^{t + to;T)=T(^]}{t;T-to). 

In particular, by setting t = and replacing the arbitrary time displacement to by t, one finds that 
the response of the medium depends only of the time difi^erence t — t, 

T^]}{t;T) = T(2iO;T-t) 
= R^^l{r-t) 

where we defined the linear polarization response function J?|ta(r — t), being a rank- two tensor 
depending only on the time difference t — t. 

To summarize, the linear contribution to the electric polarization density is given in terms of 
the linear polarization response function as 



/oo 
i?W(f-r)i^„(r,r)dT, 
-OO 



-OO 



/oo 
R(^l{T')E„{v,t-T')dT', 
-OO 



(8) 



and the causality condition for the linear polarization response function requires that 

and in addition, since the relation (8) is to hold for arbitrary time evolution of the electrical field, 
and since the polarization density P^^^(r, and the electric field Ea\r,t) both are real- valued 
quantities, we also require the polarization response function {T—t) to be a real- valued function. 
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Quadratic polarization response function 

The second-order, quadratic polarization density can in similar to the linear one be phenomeno- 
logically be written as the sum of all infinitesimal previous contributions in time. In this case, 
we must though include the possibility that not all c;ontributions origin in the same time scale, 
and hence the proper formulation of the second order polarization density is as a two-dimensional 
integral, 

/oo />oo 
/ T^^lp{t;n,T2)E^{v,Tr)E0{v,T2)dndT2. (9) 
-OO J — oo 

(2) 

The tensor TI[J^{Uti,T2) uniquely determines the quadratic, second-order polarization of the 

medium. However, the tensor T^^y,p{^^i^^2) itself unique. To see this, we may express 

the tensor as a sum of a symmetric part and an antisymmetric part, 



^;i;y*;n,T2)-T^^^i(t;r2,Ti) 



symmetric antisymmetric 



As this form of the response function is inserted into the original expression (9) for the second 
order polarization density, we immediately find that it is left invariant under the interchange of the 
dummy variables (a, n) and (since the integration is performed from minus to plus infinity 

in time). In particular, it from this follows that 



/OO nOO 
/ [4«V*; ^1' ^2) + Af^pit; ri, T2)]i?„(r, T^)Ep{v, T2) dn dr^ 
-00 J —00 



-00 ^ —00 
fOO /"OO 



/oo /"OO 
/ 'S'ylaV*; n, T2)£;„(r, Tx)Ep{v, T2) dTx dT2 
-00 J — 00 



roc poo 



+ 



/oo poo 
/ ^!S/3(*; n,T2)Eair, Ti)Ep{r, T2) dn dT2 
-00 J —00 



/OO poo 
/ '5^^ja(*; '^2, Ti)Ea{r, Ti)E(3{r, T2) dvi dT2 
-00 J — 00 

/oo pOO 
/ ^M/3a(*; 'r2,Tx)E„{Y, Ti)E0{t, T2) dTi dT2 
-00 J —00 

/oo pOO 
/ ^l^fiai*' '^2, n)Ef}{r, T2)Ea{r, n) dn dT2 
-00 J — 00 



=Il 

foo poo 



£0 



/oo poo 
/ ^Sa(*; ^2, n)E0{r, T2)E^{r, n) dn dT2 
-00 J — 00 



i. e. the symmetric part satisfy the trivial identity Ii — Ii, while the antisymmetric part satisfy 
I2 = —I2, i- e- the antisymmetric part does not contribute to the polarization density, and that 
we may set the antiymmetric part to zero, without imposing any constraint on the validity of the 
theory. 

(2^ 

The time response T^a/si^'^'^''-^'^'^) then unique and now chosen to be symmetric, 

TZ(f^n,T,) = Tlfl{t;r2,n). 
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This procedure of symmctrization may seem like an all theoretical contruction, somewhat out of 
focus of what lies ahead, but in fact it turns out to be an extremely useful property that we 
later on will exploit extensively. The previously described symmetric property will later on, for 
susceptibility tensors in the frequency domain, be denoted as the intrinsic permutation symmetry, 
a general property that will hold irregardless of whether the nonlinear interaction under analysis 
is highly resonant or far from resonance. 

By again applying the arguments of time invariance, as previously for the linear response 
function, we find that 

+ *0; n, T2) = T^'^^i; Ti - to, T2 - to) 

for all t, Ti, and T2- Hence, by setting t = and then replacing th arbitrary time to by t, again 
as previously done for the linear case, one finds that T^'^\t;Ti,T2) depends only on the two time 
differences t — Ti and t — T2- To make this fact explicit, we may hence write the response function 
as 

^S(*;n,T2) = ii£^(i-ri,i-r2), 
giving the canonical form of the quadratic polarization density as 



/■ -OC- f- OO 
/ 4a/3(i -rut- T2)E„{V, n)Ep{v, T2) dn dT2 
-OO — OO 



-OO — OO 



/OO nOO 
/ R^ifUr[,4)E^{r, t - ri)E0{r, t - 4) dri dr^. 
-OO J —CO 



(10) 



(2') 

The tensor -^^^^^('''i, T2) is called the quadratic electric polarization response of the medium, and in 
similar with the linear response function, arguments of causality require the response function to 
be zero whenever n and/or T2 is negative. Similarly, the reality condition on Ea{r,t) and Pa(r,t) 
requires that ^^Ita/jlTi, T2) is a real-valued function. 

Higher order polcirization response functions 

The nth order polarization density can in similar to the linear (n = 1 and quadratic (n = 2) ones 
be written as 



/OO POO 
■■■ / R^;2,-aAt-Tl,...,t-Tn)E„^ir,n)---E„Jr,Tn)dTi ■■■ dTn 
-00 •/ —00 



-00 
^00 



/OO roc 
■■ -R^»ai-c«„ {t[,..., T^)Ea^ (r, t - t[) ■ ■ ■ E^^ (r, t - t'J dri ■■■ dr^, 
-00 J —00 



(11) 



where the nth order response function is a tensor of rank n + 1, and a real- valued function of 
the n parameters ri, . . . , r„, vanishing whenever any < 0, /c = 1, 2, . . . , n, and with the intrin- 
sic permutation symmetry that it is left invariant under any of the n! pairwise permutations of 
(ai,Ti), . . . , (a„,r„). 
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